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^ Abstract 



We present a new theoretical approach to the kinetics of micelle formation in surfactant solutions, 
in which the various stages of aggregation are treated as constrained paths on a single free-energy 



landscape. Three stages of well-separated time scales are distinguished. The first and longest 
stage involves homogeneous nucleation of micelles, for which we derive the size of the critical 
nuclei, their concentration, and the nucleation rate. Subsequently, a much faster growth stage 

m 

takes place, which is found to be diffusion-limited for surfactant concentrations slightly above the 
critical micellar concentration (cmc), and either diffusion-limited or kinetically limited for higher 
concentrations. The time evolution of the growth is derived for both cases. At the end of the 

00 

growth stage the micelle size may be either larger or smaller than its equilibrium value, depending 



on concentration. A final stage of equilibration follows, during which the micelles relax to their 
equilibrium size through fission or fusion. Both cases of fixed surfactant concentration (closed 
system) and contact with a reservoir of surfactant monomers (open system) are addressed and 
found to exhibit very different kinetics. In particular, we find that micelle formation in an open 
system should be kinetically suppressed over macroscopic times and involve two stages of micelle 
nucleation rather than one. 
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I. INTRODUCTION 



Self-assembly of amphiphilic molecules into mesoscopic (micelles) in solution is a com- 
mon and thoroughly investigated phenomenon p|. Dynamic aspects of this process have 
been extensively studied as well [IJ. The techniques applied in such experiments and the 
interpretation of their results have used the framework of reaction kinetics, where each ag- 
gregate size is treated as a distinct chemical species, and changes in size and population 
- as chemical reactions (ref Q, chapter 3). Two well-separated time scales are identified in 
experiments [3|. The shorter of the two, denoted T\ (typically ~ 10 _5 -10~ 4 s), corresponds 
to the exchange of a single molecule between a micelle and the monomeric solution; during 
this time scale the number of micelles remains essentially fixed. The second time scale, r 2 
(which widely varies and may be, e.g., about ~ 10 -2 s) is associated with overcoming the 
barrier to the formation or disintegration of an entire micelle. The total activation time for 
such a process is mr 2 , where m is the number of molecules in a micelle. During this longer 
time scale the number of micelles changes. 

The first and still prevalent theory of micellar kinetics by Aniansson and Wall [4] is based 
on these observations. It casts the micellization process in the form of reaction kinetics with 
two well separated time scales, whereby micelles form and disintegrate through a series of 
single monomer-exchange reactions. While various extensions to the Aniansson- Wall the- 
ory have been presented over the years [5Ml3l|. only a few alternative approaches have been 
suggested. In ref the interesting possibility that micellization may behave as a bistable 
autocatalytic reaction was explored. An idealized (one-dimensional) nucleation model for 
linear aggregates was suggested in ref [H. An important alternative approach to study 
micellization kinetics has been the use of computer simulations [l6l-[23l|. In the case of mi- 
cellization of amphiphilic block copolymers more progress has been achieved (ref [2, chapter 
4; refs H3-H2) • The kinetics of such polymeric micelles, however, usually depends on qual- 
itatively different effects — in particular, the high entropy barrier for polymer penetration 
into a micelle. 

In the current work we present a new approach to the kinetics of surfactant micellization, 
which is based on a free-energy formalism. A similar strategy was previously applied to 
the kinetics of surfactant adsorption at interfaces [33l . |34| . This approach has two main 
advantages. The first is that it provides a more unified description of the kinetics — rather 
than considering different stages as separate processes ("reactions"), they can all be cast 
as constrained pathways on a single free-energy landscape. Considering different processes 
on the same footing allows, for example, easier identification of rate-limiting stages such as 
diffusion-limited or kinetically limited ones [34| . The second advantage of such a formalism is 
that it can be relatively easily extended to more complex situations, such as ionic solutions 
or surfactant mixtures [35j ]. On the other hand, the shortcoming of the model is that it 
is phenomenological, following coarse-grained thermodynamic variables rather than those 
characterizing single molecules and aggregates. It is probably not appropriate for large 
polymeric micelles, where intra-chain degrees of freedom play an important role and a more 
detailed description of molecules and aggregates is required [24], |27|, [28| . We shall focus here, 
therefore, on the micellization of short-chain surfactants. 

Another consequence of the coarse-grained modeling is that the derivation for the kinetics 
of micellization bears similarities to the kinetics of first-order phase transitions — an anal- 
ogy that was previously invoked [l5|, |27| . However, unlike macroscopic phase separation, 



micellization is restricted to finite-size aggregates, resulting, for example, in growth laws 
that are not scale-free. 

In the next section we present the free-energy formalism and its implications for the 
process of micelle formation. As in previous theories we subsequently separate the kinetics 
into stages of disparate time scales, during each of which a different set of constraints is 
imposed. We discuss separately the kinetics of closed and open systems. A closed system 
contains a fixed number of surfactant molecules. In an open system the surfactant solution 
is in contact with a large reservoir, which is at thermodynamic equilibrium. Whereas in 
equilibrium this distinction is usually immaterial, the kinetics of the two cases are found to 
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be strikingly different. While reading through the various stages of micellization it may be 
helpful to refer to the two schematic diagrams provided at the end of the article ({TBI and [TH 
for closed and open systems, respectively). The first stage that we address is the nucleation 
of micelles. Subsequently, we describe the growth of the micellar nuclei as they absorb 
additional monomers from the surrounding solution. Both options of kinetically limited and 
diffusion-limited growth are studied. In addition, the possible role of long-distance diffusive 
transport is examined. We then consider the final relaxation toward equilibrium. Finally, 
we summarize the conclusions and discuss the experimental implications of our analysis, as 
well as its limits of validity. 



II. MODEL 

The model is based on a simple free-energy functional, which has been recently introduced 
to study metastability issues of micellization [36|. Apart from the temperature T, the 
free energy depends on three thermodynamic degrees of freedom, which we take to be the 
total volume fraction of surfactant in the solution, $, the volume fraction of surfactant 
monomers, $1, and the number of molecules in a micelle (aggregation number), m. Despite 
the simplified two-state (monomer-aggregate) description, polydispersity can be accounted 
for as fluctuations of the variable m [33] . (This, however, restricts the validity of the model 
to compact micelles whose size distribution is narrow p].) All energies hereafter are given 
in units of the thermal energy, k-^T. 

The free energy has contributions from the entropy of mixing and from the interactions 
among surfactant molecules. The former is obtained from a coarse-grained lattice scheme 
(Flory-Huggins model), in which a water molecule occupies a single lattice cell of volume a 3 , 
and a surfactant molecule is larger and occupies n such cells. The interactions in the solution 
are represented by a single phenomenological function, u(m), which is assumed to capture 
all the molecular contributions to the free energy of transferring a surfactant molecule from 
the solution into an aggregate of size m. The resulting Helmholtz free energy density (per 
lattice site) is [36| 

F($,$i,m) = — ln$i + — [In $ m - mu{m)\ + (1 - $)ln(l - $), (1) 

where $ m = $ — $i is the volume fraction of micelles, and 1 — $ is the volume fraction of wa- 
ter. At equilibrium the solution is spatially uniform and characterized by those single mean 
values of the variables, which minimize the free energy under the appropriate constraints. 
For a closed system F is minimized with respect to $i and m for a given $. For an open 
system one should minimize F — /i$ with respect to $, $i, and m for a given surfactant 
chemical potential /i. Out of equilibrium the values of variables, such as $, $i, and m, 
may be position-dependent, and the total free energy is given by spatial integration of the 
local free-energy density. (We neglect here surface-tension (gradient) terms associated with 
boundaries between such spatial domains.) 

The specific choice of the interaction function u(m) is not crucial so long as it has a 
maximum at a finite m to ensure the stability of finite-size micelles. To provide numerical 
examples, and following previous works [3g, |38[, we use a simple three-parameter function, 

u(m) = uq — am^ 1 ^ — nm 2 ^ 3 . (2) 

The first term in eq [2] represents a micelle-size- independent free-energy gain in increasing m, 
the second — a surface energy penalty, and the third is responsible for stabilizing a finite- 
size aggregate. (For a more detailed discussion of these terms and the restricted ranges of 
relevant values for Uq, a, and k, see ref l36l.) 

Despite its simplicity eq [I] defines a rather rich free-energy landscape over a three- 
dimensional space of macro-states, parametrized by the axes ($, $i,m). We are going to 
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treat the kinetics of micellization as time- dependent paths along this landscape, and it is 
beneficial, therefore, to first recall its key features [36f| , demonstrated in[TJ For any given <3> 
and along the $1 axis, F always has a single minimum at $1 = $^(m, $) for all values of 
m. This value of monomer volume fraction as a function of aggregation number and total 
volume fraction is found by solving the equation 

= $*(m, $) : ^m e mu(m)+m-l = $ _ ^ (3) 

Along the m axis, however, F becomes nonconvex when (f) exceeds a certain volume fraction, 
tpi, with two minima at m — 1 and m = m*($i, $), and a maximum in between at m — 
^nuc($i, (See[T] A and B.) The extrema satisfy the equation 

m = m* , m nuc : m 2 = — ln($ — '(m), (4) 

where w' = du/dm. Combining eqs [3] and HI we can find m and $i at the extrema for a 
given $ according to 

m = m*,m nuc : m 2 = - ln[$ - e-" (m) ~ m "' (m) - 1+1/m ]/u '(m), (5) 

$* = e -«H-w'H-l+l/m^ ^ 

Above a larger volume fraction, ^ > fi (IIP), the micellar state with $ > </?2, Tn = m* 
and $i = $J(m*,$), though still metastable, may become appreciably occupied, giving 
rise to premicellar aggregates [36]. Above yet another volume fraction, tp 3 > (p% (pD), the 
micellar state for $ > <^ 3 becomes the global minimum of F . It is this point, analogous to 
the binodal line in phase separation, which corresponds to the commonly defined cmc [36], 
i.e., cmc = c 3 = (p 3 /(na 3 ). We shall focus in the current work on the ordinary micellization 
region, $ > tp 3 = ip crac , where micelles are stable at equilibrium. It should be kept in mind, 
however, that in this region the monomeric and micellar states are separated by a free-energy 
barrier in the form of a saddle point of F, F nuc ($) = F[$, <3>^(m nuc , $), m nuc ]. The barrier 
may be high, leading to the measurement of an apparent cmc, which is higher than the 
equilibrium one, tp 3 = feme [El- Finally, above a certain higher volume fraction, y? 4 > tp 3 , 
the barrier disappears and the micellar state for $ > y9 4 remains the sole minimum of F, as 
seen in [T]E. (This work is restricted to the isotropic micellar phase of surfactant solutions; 
at higher concentrations other phases and meso-phases appear The point $ = if4 IS 
the analogue of the spinodal line in macroscopic phase separation. As already mentioned 
in the Introduction, despite the analogy with phase separation it should be borne in mind 
that micellization is essentially different in that it involves finite-size aggregates and smooth 
crossovers rather than macroscopic phases and sharp transitions. 

The initial and final states of the micellization kinetics are defined as follows. At t — 
the system is in the monomeric state, ($i = $,m = 1), whereas its equilibrium state is the 
micellar one. In a closed system this is done by setting the surfactant volume fraction above 
the cmc, $ > f 3 (using, for example, the temperature-jump or stopped-flow techniques 
[2]). In an open system the initial condition corresponds to opening a diffusive contact with 
a bulk reservoir, whose surfactant volume fraction <3>b is above the cmc, <3>b > y?3- The 
reservoir is assumed to have already reached the equilibrium micellar state. At t — > oo the 
system reaches the global minimum of the free energy — [$, $J(m*, $), m*($)] in the closed 
case and [$b, 3>i(m*, $b) ; ^(^b)] m the open one. In what follows we consider the kinetic 
pathway that the system takes between these initial and final states. Assuming separation 
of time scales, we shall divide the temporal path into separate stages. Note that the various 
time scales are derived from the free energy functional and a single molecular time, r , thus 
enabling comparison of different stages and processes. 

Throughout the following sections we demonstrate the results using a sin gle exemplary 
surfactant, whose parameters are listed in[H This allows comparison with refs 13a. I37I. where 
the behavior of the same exemplary surfactant for $ < ^3 = <p cmc was presented. 
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FIG. 1: Cuts of the free-energy landscape (per surfactant molecule, in units of k^T) as a function 
of aggregation number along the §\(m) line for the surfactant parameters of U and increasing 
surfactant volume fraction, (A) $ = 5 x 1(T 4 < ipx; (B) (pi < $ = 7 X 10~ 4 < (p 2 ; (C) 
ip 2 < $ = 1.5 x 1CT 3 < <^3 = v? cmc ; (D) 993 < $ = 2.5 x 1CT 3 < <^ 4 ; and (E) $ = 0.11 > </? 4 . 

[2] shows two cuts through the free-energy landscape as a function of aggregation number 
for the exemplary surfactant in a closed system at total surfactant volume fraction slightly 
larger than (p% = y2 cmc . Along the first cut (solid line) the monomer volume fraction is 
assumed to be at quasi-equilibrium, $1 = $*(m). Thus, the minimum of this curve corre- 
sponds to the global minimum — the equilibrium aggregation number. Along the other cut 
(dashed curve), which is relevant to the next two sections, we constrain the concentration 
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TABLE I: Parameters of the exemplary surfactant and the resulting boundaries of the micellar 
region 



n u Q a k tp 3 = ip cmc ip 4 



13 10 11 0.08 2.03 x 10~ 3 0.106 
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FIG. 2: Cuts of the free energy landscape (relative to the monomeric state, per micelle, in units of 
k&T) as a function of aggregation number for the surfactant parameters of [J and $ = l.lc/23. The 
two curves correspond to two different constraints: relaxation of the monomer volume fraction for 
the given aggregation number (solid), and fixed concentration of micelles (dashed). Indicated by 
arrows are the sizes of the critical nucleus (m nuc ), the intermediate aggregate at the end of the 
growth stage (m), and the equilibrium micelle (rn*). A closed system is assumed. V\ is the volume 
of solution per micelle at equilibrium. 



of micelles to remain at its nucleation value. 



III. MICELLAR NUCLEATION 



A. Closed System 

Let us set the total volume fraction at t = to some value, $ > (f 3 = ip cmc , and 
assume that value (apart from a short initial period of homogenization which is ignored) to 
remain fixed and uniform throughout the micellization process. The first stage to consider 
is the ascent of the free energy from the initial metastable state, ($, $1 = $,m = 1), to 
the saddle point [$, $1 = §l(m nuc ),m = m nuc ] - i.e., the formation of the critical nuclei. 
This activated process is assumed to be much slower than diffusion. Hence, $1 can be 
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taken during this stage as spatially uniform and equal to the value that minimizes the free 
energy for the given $ and m(t). Thus, as m(t) increases from 1 to the critical-nucleus 
size m QUC , the system proceeds along the path that satisfies the constraints $ = const and 
$ x = $*[m(t),$]. 

A similarly constrained path was studied in detail in ref |37] to obtain the lifetime of 
metastable micelles in the region <p2 < $ < ^3 using Kramers' theory. Such a rigorous 
calculation, unfortunately, cannot be repeated here, since the metastable monomeric state 
is actually not a local minimum of F but just the edge, at m = 1, of the range of allowed 
aggregation numbers. (See[Tp.) Nevertheless, as demonstrated in refH3, the nucleation time 
(dissociation time in ref [37T) and its concentration dependence are primarily determined by 
the height of the free-energy barrier. 

The free-energy barrier corresponds to the nucleation of a single micelle. Our model, 
however, considers macrostates of a solution containing many micelles and monomers. To 
switch between these two descriptions we introduce a subsystem volume, Vi, which contains 
(on average) a single nucleus. The volume fraction of critical nuclei, their concentration, 
and the volume per nucleus are readily given for closed systems by 

$ nuc ($) = $-<^[m nuc ($),$] 
($) = $ nuc ($)/[na 3 m nuc ($)] 

where m nuc ($) and <&*($) are given by eqs El and [6] Since $ nuc is very small, V\ is much 
larger than the molecular volume, and our coarse-grained approach is indeed applicable. 
Note the distinction between the nuclei concentration c nuc and their volume fraction $ nuc . 
Since the micelle size m is a variable, constraining c nuc does not imply a fixed $ n uc- This will 
be important in the next sections, when we impose a constraint on the number of nuclei. 
The nucleation barrier and nucleation time scale are given for closed systems by 

Tnuc($) * T e AF ™W, (8) 

where Tq is a molecular time scale, and F\ is the free energy of the monomeric state. It should 
be mentioned that our formalism artificially distinguishes between monomers and aggregates 



of size m = 1 . As in the previous works |36|, [37J , this artifact has an insignificant effect on the 



results. We calculate here the free energy of the m = 1 state as i r i($) = F[<&, $*(1, $), 1]. 

Various features of the nucleation stage can be calculated from eqs dHEJ as demonstrated 
in Figures EH5J The concentration of critical nuclei ([3A) sharply increases with surfactant 
volume fraction as $ is increased above <£> 3 = ^) cmc . The size of the critical nucleus (TJA) de- 
creases with $ until it practically vanishes as $ approaches (p^. The height of the nucleation 
barrier (jHJ) decreases as well with $, leading to a sharp decrease in the nucleation time scale 
([5] inset). To get an estimate of the actual nucleation time scales we may take r ~ 10 -8 
s, which is the time it takes a molecule with a diffusion coefficient ~ 10 -6 cm 2 s" 1 to be 
displaced by ~ 1 nm. For the example presented in El T nuc is extremely large close to <^ 3 but 
drops to ~ 1 s for $ ~ 2<^3. 



B. Open System 

When the system is placed in contact with a large reservoir of volume fraction $b, $b > 
^3 = V^cmc; monomers will first diffuse in, until the monomeric concentrations are balanced. 
We shall assume that micellar diffusion from the reservoir is either blocked or very slow. (If 
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FIG. 3: Concentration of critical nuclei (normalized by the molecular volume) as a function of 
surfactant volume fraction in the range between (^3 = ip cmc - 2x 1CF 3 and ip4 for closed (A) and 

□ 




FIG. 4: Critical-nucleus size as a function of surfactant volume fraction in the range between 
(£>3 = (feme — 2 x 10~ 3 and ^4 for closed (A) and open (B) systems. Parameters are given in|TJ 

it is not, micellization in the system will be dominated by simple transport of micelles from 
the reservoir.) Thus, the starting point for the nucleation stage in this case is different from 
that of a closed system — it is still a monomeric state, yet with a lower volume fraction, 
$ = $1 where $1 = &\ < $b- Nucleation is again assumed much slower than monomer 
diffusion. Hence, the monomer volume fraction remains fixed at $1 = At the same time 
it should minimize F for the given m(t), which in turn determines the value of the third 
state variable, $. As the nuclei grow, the total volume fraction increases, and the system 
proceeds along the path that satisfies the constraints $1 = &\ and $*[m(t), $] = $>\. 

The nucleation path ends at the state of critical nuclei, which is also different from the 
closed-system saddle point, because the total volume fraction has not reached the bulk value, 
$ < $b- This state is calculated using the following procedure for the open case. First, 
we calculate the monomer volume fraction in the reservoir according to the equilibrium 
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FIG. 5: Nucleation barrier A.F nuc (in units of k^T) as a function of surfactant volume fraction for 
a closed system. The inset shows the corresponding nucleation time (in units of the molecular time 
To). Parameters are given in [J 

condition, 

$5 , ($ b ) = $ 1 K($ b ),$ b ]. (9) 

Second, we equate this monomeric volume fraction with the one in our open system at the 
saddle point, 

$*K uc ($),$] = $t($ b ), (io) 

thus determining (via eqs |5] and |6]) the total volume fraction in the system, $, and the 
critical nucleus, m nuc , as functions of $b- Third, we use these results to calculate $ n uc, c nuc , 
and V\ as functions of $b) 

$nuc($b) = $ - $1 

Cnuc($b) = $nuc/(na 3 m 



nuc ) 

3„ 



_ j na m 



Finally, the nucleation barrier and time scale are given for the open system by 

AF nuc ($ b ) = ^[F^^mJ-f!^)] 

W($b) ^ r e AF -(* b ). (12) 

From eqs dHS] and l9"HT2l one can calculate the various parameters of the nucleation stage 
for an open system. Examples are shown in Figures [3)3 andH]B, revealing striking differences 
from the case of a closed system. The explanation is straightforward — the system is assumed 
to be in contact with the reservoir only through its monomeric concentration (so-called inter- 
micellar concentration), <&\, which hardly changes as $b is increased above the cmc. Hence, 
during this initial stage <3>i remains low regardless of the value of $b- Consequently, the 
critical nuclei remain relatively rare and large, almost independent of concentration (Figures 
[3)3 and Hf3). Moreover, since $i does not reach values above y? 3 = y3 cmc , we get very high 
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nucleation barriers, resulting in an unphysical nucleation time for the open system. Thus, 
homogeneous nucleation in an open system, which does not have micellar transport from 
the reservoir, is strongly hindered. In the following discussion of open systems it is assumed 
that, despite this strong kinetic limitation, nuclei were somehow caused to form. 



IV. MICELLAR GROWTH 

The nucleation stage addressed in the preceding section ends when the critical nuclei 
have formed. On the free-energy landscape the system has reached the saddle point of F. 
Subsequently, a stage of faster growth takes place. The free energy of the system decreases 
while the nuclei absorb additional monomers from the surrounding solution and m increases. 

The growth is assumed to be much faster than the nucleation of new micelles or fusion 
and fission of existing ones. Hence, the concentration of micelles, c m = ($ — $i)/(na 3 m) 
remains fixed at c m = c nuc . Consequently, the available volume per aggregate, V%, remains 
unchanged as well. We shall assume that the growth is also faster than the diffusive transport 
among the micelles (for closed and open systems) and with the reservoir (open system). The 
increase in m, therefore, comes solely at the expense of a decrease in the concentration of the 
surrounding monomers, while the total surfactant volume fraction is conserved. Thus, we 
describe the growth kinetics as a constrained path, [3>i(t), Tn(t)), such that c m = c nuc = const 
and $ = const. 

Although diffusive transport into or out of the subsystem (of volume V\) is assumed neg- 
ligible during this stage, it is a priori unclear whether the growth process itself, within V\, 
should be kinetically limited or diffusion-limited. We shall therefore examine both possibil- 
ities below. The constraints and the equations derived in this section apply to closed and 
open systems alike, yet the values substituted for $ and c nuc differ substantially. While for 
a closed system $ is the experimentally controlled surfactant volume fraction, for an open 
system <3> gets the lower and weakly changing values determined from $b i n the nucleation 
stage according to eq[]IJJ The concentration of nuclei is also much lower in the open-system 
case (cf. |3]). Consequently, the quantitative predictions for the two cases are quite different. 

The aforementioned constraints imply that the average monomer volume fraction de- 
creases linearly with the aggregation number, m(t), 

= $ - na 3 c nuc m(t). (13) 

We are left with one independent variable, m(t), whose change in time could be either 
kinetically controlled or diffusion-controlled. Yet, before studying the detailed evolution, let 
us examine its final state, which is common to both limits. 

The final state of the growth stage, denoted (5> 1; m), is given by the minimum of F along 
the constrained path, (dF/dm) | =c $ =comrt = 0. This yields 

m = rh : ln[$i(m)] + u(m) + mu'(m) + 1 — 1/m — 0, (14) 

where $i(m) is given by eq [131 an d, once rh is calculated, $i = $i(m). The resulting 
aggregation numbers and their dependence on the controlled surfactant volume fraction are 
presented in [6j Note that the intermediate aggregation number at the end of the current 
stage is not equal to the equilibrium micellar size, since it corresponds to a minimum of F 
along the constrained path rather than its global minimum. Unlike the equilibrium size, 
m*, which is bound by thermodynamic stability to increase with surfactant volume fraction 
(dotted lines in [6]), the intermediate size rh can have a richer behavior. Examined over a 
wider range of $, rh is found to be nonmonotonous, having a maximum at $ < (p^ = y? cmc . 
Hence, for the closed system it decreases with $ (jSK), whereas for the open system, which 
remains dilute throughout this stage, it increases with $ (and, therefore, with <&t,; HUB). In 
the closed system the growth overshoots the equilibrium size for $ > y? 3 and undershoots it 
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FIG. 6: Intermediate micelle size at the end of the growth stage, 111 . duS £L function of surfactant 
volume fraction in the range between (^3 = (^ cmc ~ 2 x 1CP 3 and tp^ for a closed (A) and open (B) 
systems. The inset in panel A focuses on volume fractions slightly above 993. Dotted lines show 
tl e □ 




FIG. 7: Intermediate monomer volume fraction at the end of the growth stage as a function of 
surfactant volume fraction in the range between ^3 = i/3 cmc ~ 2 x 10~ 3 and ^4 for a closed (A) 
and open (B) systems. The inset in panel A focuses on volume fractions slightly above ip%. Dotted 
lines show for comparison the equilibrium monomer volume fraction, Parameters are given in 

m 

at higher values. Whether rh is larger or smaller than m* is in accord with the question of 
whether c nuc is, respectively, smaller or larger than the equilibrium concentration of micelles. 
(We shall return to this point when we deal with the final relaxation.) In the open system rh 
is very close to, and slightly smaller than, m*. Similar observations can be made concerning 
the intermediate monomer volume fraction, $1, as demonstrated in [71 

We now turn to the evolution of the micellar size. We shall first assume, in the first 
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FIG. 8: Temporal increase in micellar size assuming kinetically limited growth in closed (solid line) 
and open (dashed line) systems. The curves are obtained from numerical solution of eq[15]for the 
parameters given in [J a = 1, and = 2.23 x 10~ 3 = for the closed system, while for the 

open system the same value is taken for <l>b- 

subsection below, that it is kinetically limited. We will subsequently check in the second 
subsection whether such a description is consistent with the rate of monomer diffusion and 
consider the alternative of a diffusive growth. 



A. Kinetically Limited Growth 

In the case of kinetically limited growth the diffusive transport of molecules to the ag- 
gregate is assumed sufficiently fast so as not to limit the growth. The volume fraction of 
monomers, $1, satisfies eq [13] while being uniform across the subsystem volume V\. The 
increase of m with time is taken as proportional to the relevant thermodynamic driving force 
(i.e., the slope of F along the constrained path), 

ot 

= — {ln[$i(m)] + u{m) + mu (m) + 1 — 1/m} , (15) 

Cm — C nuc 

$=const 

where a is an unknown dimensionless prefactor of order unity, and $i(m) is given by eq 
[T3l Equation [151 supplemented by a proper initial condition for m(t = 0), forms a simple 
initial-value problem for the temporal increase in micelle size, and is solved numerically. 
Since the initial state of this stage is a stationary (saddle) point of F, we cannot begin 
with the strict initial condition, m(0) = m nuc , but have to perturb it to start the evolution. 
An example for a numerical solution of eq[TS], where we have taken m(0) = m nnc + 1 and 
$ = l.lip 3 = l.lip cmc , is shown in[HJ The time scale of the growth, denoted r k , is found to 
be about two orders of magnitude larger than the molecular time r (i.e., of order 10 -6 s in 
this example). 

To get an expression for the kinetic time scale we examine the asymptotic behavior of eq 



dm a V\ 5F 
dt r a 3 5m 
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FIG. 9: Time scale of kinetically limited growth as a function of surfactant volume fraction in the 
range between 993 = ^3 cmc ~2x 10~ 3 and ^4 for closed (A) and open (B) systems. The insets show 
the increase of t^T 1 (growth rate) with <I> (in A) or <I>b (in B). Parameters are given in [J and we 
have set a = 1 in eq [16j 



IT5las m approaches m, obtaining 



\m(t) 



m ~ e 



-t/Tk 



Tu = — 



a 

TO 



$ - $1 



2u '{m) — mu "{m) — 1/ m 2 



(16) 



m=m, &i=&i 



The results for Tk in terms of the molecular time tq are shown in |9j For the closed system, 
over one decade of surfactant volume fraction, decreases from ~ 10 2 r to ~ r . (Values 
below t , evidently, should not be regarded as physical.) The inset shows that the growth 
rate for the closed system increases roughly linearly with $. For the open system the time 
scale is also about two orders of magnitude larger than r , yet its dependence on $ b is much 
weaker for the reasons described in the Nucleation section. 



B. Diffusion-Limited Growth 



In the preceding subsection we have assumed that the surrounding solution can supply 
the amount of monomers required for micellar growth within the time scale t^. Let us check 
whether this assumption is consistent with the rate of diffusive transport from the solution 
into the aggregate. The thickness of the diffusion layer around the aggregate, £4 (assumed to 
be much larger than the aggregate radius), satisfies the equation Am = (47r/3)£^ci, where 
Am = fh — m nuc is the number of monomers to be transported, and c\ = $i/(na 3 ) the 
monomer concentration. The diffusion time scale is then ~ D being the diffusion 

coefficient of a monomer. Using the definition to ~ (na) 2 /D, we obtain 

r d /r ~ [3Am/(4v™ 2 )] 2/3 $^ 2/3 ~ (0.1-1)<^ 2/3 , (17) 

where in the last relation we have assumed n ~ 10 and Am ~ 50. For our typical example 
of $1 ~ 10~ 3 (cf. [7]), we get r d ~ (10-10 2 )r , i.e., comparable to r k . Thus, the situation 
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concerning the limiting process for micelle growth is not clearcut, and both processes may 
be relevant in general. 

To treat the diffusion-limited growth in more detail we employ the following approxima- 
tions. First, we neglect the increase in the aggregate radius, R, and take it as constant. 
Although this description is evidently inaccurate, it crucially allows us to avoid the com- 
plicated treatment of a moving boundary. Since the growth does not begin from a single 
monomer but from a critical nucleus of finite size m nuc , we do not expect the approximation 
of constant R to qualitatively affect the results. Second, the diffusion layer is assumed much 

1/3 

smaller than the subsystem, -C V 1 , thus allowing us to consider the latter as infinite, 
and the monomer volume fraction far from the micelle as given by eq [13j Third, we neglect 
desorption of monomers from the micelle to the solution during the growth. This is justified 
in view of the strong driving force (large slope of F) for growth above the critical-nucleus 
size. 

We assume a radial volume- fraction profile of monomers, $i(r > R,t), which follows the 
diffusion equation, 

<9$! n l 9 { 2 9 ®1 



dt r 2 dr \ dr ) ^ ^ 

The growth of a micelle is determined by the diffusive flux of monomers from the solution, 



dm jj^ R2 



(19) 

r=R 



dt na 3 dr 

The boundary condition far from the micelle is given according to eq [13] by 



$i (r ->• oo,t) =<& — na 3 c nuc m(t). (20) 

For the problem to be well posed, eqs [1814201 should be supplemented by appropriate initial 
conditions for $i(r, 0) and m(0), as well as a local "adsorption isotherm" at the aggregate 
surface, relating &i(R, t) and m[t). The latter lies beyond the scope of our coarse-grained 
description. At any rate, we are interested primarily in the qualitative asymptotics of the 
diffusive transport from the solution into the aggregate, for which these details are not 
crucial. The asymptotic behavior as the final state of the growth stage is approached is 
worked out in the Appendix, yielding 



^(R, oo) - [1 - (r d /tp] , rd = a2( ^ )2/3 $r 2/3 - (21) 

Thus, unlike the exponential relaxation of a kinetically limited process (eq [T6|) . the diffusive 
relaxation is characterized, as usual, by a slow power law. Upon substituting r ~ {na) 2 / D 
in eqEDthe general form of Td, derived earlier from heuristic arguments (eq lTTj) . is confirmed. 

[TOl shows the dependence of Td on the controlled surfactant volume fraction according to eq 
[2T1 where we have taken r = (na) 2 / D. The cases of closed and open systems are again found 
to behave qualitatively differently, strongly decreasing with $ in the former and weakly 
increasing with <3>b i n the latter. This is a consequence of the different dependencies of rh on 
concentration, commented on earlier (cf. [6]). In an open system rh increases with $b P^B)] 
and, since the more molecules are transported the longer the diffusive process takes (i.e., 
increases with Am in eq l2Tl) . we get an increase of t& with $b [TOT B)]. Comparison of Figures 
191 and [TUl confirms our earlier assessment, that and Td are comparable in general, and both 
growth mechanisms may be relevant. Only for a closed system at concentrations slightly 
above the cmc do we get for our representative example Td 3> T\, i.e., strictly diffusion- 
limited growth. (Note that and Td are associated with very different time dependencies 
- an exponential law vs. a power law — and are defined only up to a numerical prefactor. 
Hence, they should be compared with respect to the order of magnitude only.) 
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FIG. 10: Time scale of diffusion-limited growth as a function of surfactant volume fraction in the 
range between ip% = ^ cmc ~ 2 x 10~ 3 and 934 for closed (A) and open (B) systems. Parameters are 
given in|Il 

C. Role of Bulk Diffusion 

In the preceding subsection we have considered the local diffusive transport that takes 
place around individual micelles, feeding them with monomers. In the case of an open system 
there should also be slower, long-distance diffusion of monomers from the bulk reservoir. In 
principle this should have been the next stage to consider. However, we find that the 
monomer volume fraction at the end of the growth stage, $1, is invariably very close to 
the equilibrium (bulk) value, $*. (See [33-) This is a consequence of the small number of 
initial nuclei (IHB), whose growth consumes a small number of monomers. Thus, the driving 
force for bulk diffusion is very weak. Consistently, for the open system we find also that 
the micellar size at the end of the growth stage, m, is very close to the equilibrium size, 
m* (PB). Therefore, the bulk diffusion that does occur in an open system has a very minor 
contribution to the micellization. 



V. FINAL RELAXATION 

At the end of the growth stage monomer transport into the existing micelles has been 
exhausted, and the micelles have equilibrated with the surrounding monomers. Yet, the final 
state of this stage, ($1, fh), does not correspond to the global free-energy minimum, since up 
till now we have constrained the concentration of micelles to remain at its nucleation value 
(cf. [2]). A slower process should ensue, therefore, during which the size and/or concentration 
of micelles relax to their equilibrium values. 

In the open system the situation is a bit unusual. (Recall from the Nucleation section, 
however, that actually reaching the current stage in an open system should already involve 
overcoming unusually high barriers.) The monomer volume fraction has equilibrated with 
the bulk reservoir and reached its equilibrium value. The size of the existing individual 
micelles has equilibrated as well. What has not equilibrated yet is the total surfactant 
volume fraction — specifically, the contribution to $ from $ m , the micellar volume fraction. 
Since there is no thermodynamic driving force for either monomer transport or changes in 
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the size of the existing micelles, and we do not allow for transport of micelles from the 
reservoir, the only open pathway to final relaxation is the very slow nucleation of additional 
micelles. The newly formed micelles will take monomers from the solution, causing transport 
of additional monomers from the reservoir, until the total surfactant volume fraction reaches 
its equilibrium value, $b- 

The relaxation of the closed system is qualitatively different. Both the monomer volume 
fraction and aggregation number have not equilibrated yet and will change in time while 
keeping the total surfactant volume fraction fixed. Since there is no longer a driving force 
for directional exchange of monomers with the solution, we expect these changes to occur 
through fusion or fission of micelles. Such processes occur on the scale of an entire micelle 
and depend, therefore, on a different microscopic time, denoted r m . It is expected to be 
much larger than the molecular time tq — either because of the long diffusion time required 
for two micelles to meet before fusing (in which should be of order, say, 10 5 -10 4 

s), or due to kinetic barriers for fusion or fission. In addition, r m should depend on details 
of inter-micellar interactions. Such kinetic barriers and interactions are not accounted for 
by the current model. Two additional processes, which in principle could be considered, 
are irrelevant in this case. First, nucleation of new micelles or complete disintegration of 
existing ones might occur but will require the much longer time scale of r nuc discussed earlier. 
Second, Ostwald ripening — a common relaxation mechanism in phase separation, where 
larger domains grow at the expense of smaller ones — is not expected to take place, since 
the finite domains here (the micelles) are not unstable and the required positive feedback is 
thus lacking. 

Either fission or fusion should be dominant, depending on whether m has overshot or 
undershot, respectively, the equilibrium size m*. (See E]A..) Correspondingly, the micellar 
concentration c m will either increase or decrease with time. Over the time scale of these 
rearrangements of aggregate size and concentration we can assume that the monomer volume 
fraction is relaxed, $i(i) = $^[m(t), $]. We are left again with a single kinetic variable - 
either m(t) or c m (t). The two are related via 



c m {t) = {$ - $![m(t), $}}/[na 3 m(t)}. 
The kinetic equation for the micellar size reads 
dm (3 V\ 



(22) 



dt 



-f{m) 



f{m) = — 
dm 



<&!=$* (m) 
<J>=const 



$t'm$i - 



m c 



m 



ln($ - $*) 



(23) 



-($ - $*)«'(m) + [u(m) + 1 - l/m]$*', 



where V\ = na 3 m*/[$ — <&*(ra*)J is here the volume per micelle at equilibrium, is 
given by eqO a prime denotes d/dm, and f3 is an unknown dimensionless prefactor of order 
unity. 

Equations 131 and [231 are solved numerically to obtain m(t) and, subsequently (via eq l22l) . 
also c m (t). [TT1 shows the solutions for our exemplary surfactant and two volume fractions, 
corresponding to fission- and fusion-dominated relaxation. 

To find the relaxation time we examine the asymptotic behavior of m(t — > oo) according 
to eq[23l obtaining 



\m(t) — m* 



~ e 
1 



-t/T T 



13 V l f'(m* 



(24) 
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FIG. 11: Evolution of micellar size (solid, left ordinate) and concentration (dashed, right ordinate) 
during the final relaxation stage in a closed system. Parameters are given in HJ we have set (3 = 1 
in eqE3 and the volume fraction is $ = 2.23 x 1(T 3 = 1.1^3 (A) and 5.23 x 10 -3 = 2.58^3 (B). 
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FIG. 12: Time scale of final relaxation, r r (in units of the single-micelle time scale), as a function of 
surfactant volume fraction in the range between ip^, = ^ cmc ~ 2 x 10 -3 and ip^ for a closed system. 
Parameters are given in HI and we have set f3 = 1 in eq 1241 

where f(m) has been denned in eq[23j The dependence of r r on surfactant volume fraction 
is shown in[12j The relaxation time is found to weakly depend on remaining of the same 
order as (or slightly larger than) the single-micelle time r m throughout the concentration 
range. 
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VI. DISCUSSION 



The detailed picture, which arises from our analysis of micellization kinetics, is schemat- 
ically summarized in Figures [T3] and [HJ We have divided the process of micelle formation 
into three major stages — nucleation, growth, and final relaxation. On the one hand, this 
crude separation into stages should be conceptually valid, since we find the corresponding 
time scales to be quite well separated. In particular, the nucleation time is found to be 
macroscopic, several orders of magnitude longer than the time scales of growth and equi- 
libration. Such three stages have been resolved in a recent x-ray scattering experiment on 
block copolymer micellization [31]. They also emerged in other micellization theories [151. 
On the other hand, the discreteness of these stages, as illustrated in Figures [13] and [T4l 
should be taken too strictly. In particular, in the example treated above we find the time 
scale of growth to be only 1-2 orders of magnitude shorter than the typical time for final 
equilibration. Thus, in certain cases it may well be that these two stages should not be 
considered as distinct. 

The nucleation stage is much longer than all others and, since it is an activated process, its 
duration is exponentially sensitive to surfactant volume fraction as well as other parameters 
(J5]). The range of nucleation times that we get for our exemplary surfactant in a closed 
system (typically larger than 1 s) is in line with measured values of mr 2 — the time scale 
for formation or disintegration of entire micelles j2|. The high nucleation barriers found 
close to the equilibrium cmc ($ = ^3) imply that the measured (apparent) cmc might in 
certain cases be higher than the equilibrium value. This issue, which was raised before in 
the context of block copolymer micelles [l3| , clearly merits further study. 

The growth stage occurs on much faster time scales (e.g., 10~ 6 -10 -5 s for our example). 
These time scales are similar to those measured for t\ — the single-monomer exchange time 
at equilibrium [2|. We have found that the growth may in general be either diffusion-limited 
or kinetically limited, and that it should be diffusion-limited at concentrations close to 
the cmc. This is in accord with n being usually diffusion-limited for short-chain surfactants 
while becoming kinetically limited for longer-chain ones, which face higher kinetic barriers for 
incorporating into a micelle j2j. Our theory predicts a distinctive —3/2 power-law relaxation 
in the case of diffusion-limited growth (eq |2"TT) . This prediction should be verifiable in 
scattering experiments like the one described in ref [M], when they are applied to short-chain 
surfactants. 

The final relaxation stage in a closed system may involve either reduction in aggregate size 
(fission), accompanied by an increase in aggregate concentration, or the other way around 
(fusion). (See El) Which of these scenarios holds depends on whether the aggregate size 
attained in the preceding growth stage has overshot or undershot the equilibrium aggregation 
number. The former should hold at concentrations close to the cmc, whereas the latter occurs 
at higher concentrations. We note that in the experiment of ref [Ml the aggregates grew in size 
during their final relaxation, which is in line with the fact that the amphiphilic concentration 
in that experiment was much higher than its cmc. We note also that the somewhat surprising 
possibility of an intermediate aggregate size overshooting the equilibrium value was already 
pointed out in an earlier study [27| . An interesting consequence of our analysis is that, 
by tuning to the right surfactant concentration, one should be able to eliminate the final- 
relaxation stage altogether, thus reaching the equilibrium micellar state already at the end 
of the fast growth stage. Another relevant prediction is that the relaxation time of this final 
stage should be almost independent of surfactant concentration f[T2l . It should be stressed 
again that these predictions concerning the final relaxation stage require that the preceding 
growth stage be sufficiently fast so that the two processes could be considered separately. 
In particular, observing oversized micelles before they shed their extra molecules may be 
experimentally challenging. 

Our findings concerning the kinetics of micelle formation have a number of additional 
experimental implications. A particularly clearcut one relates to micellization in an open 
system — a solution in diffusive contact with a reservoir of monomers and micelles. We 
have found that, in cases where only monomer exchange with the reservoir is allowed while 
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FIG. 13: Schematic summary of micellization in a closed system. The states of the system are 
represented by rectangles, beside which the values of the state variables are indicated. The process 
is divided into three stages, represented by arrows. The constraints on the kinetics during each 
stage are indicated beside the arrow, (i) Slow nucleation stage (time scale r nuc ), in which critical 
nuclei of size m nuc form in a monomeric solution, (ii) Fast growth stage, in which the nuclei grow 
from m nuc to an intermediate size fh without changing their concentration. The growth may be 
either diffusion- limited (time scale t^; —3/2 power-law relaxation) or kinetically limited (time scale 
7k; exponential relaxation), (iii) Final relaxation of the size and concentration of aggregates to 
their equilibrium values through fusion or fission (time scale r r ). 



the transport of micelles is blocked, micellization should be kinetically suppressed. The 
suppression is two- fold. First, strong activation is required for the homogeneous nucleation 
of the first micelles. This stems from the low surfactant concentration maintained in the 
system due to the correspondingly low monomer concentration (sometimes referred to as 
the inter-micellar concentration) in the reservoir. Second, even after micelles do nucleate 
and grow, the final relaxation of their concentration should be hindered, since it requires 
the nucleation of additional micelles. 
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FIG. 14: Schematic summary of micellization in an open system, having monomer exchange with 
a reservoir. The states of the system are represented by rectangles, beside which the values of the 
state variables are indicated. The process is divided into four stages, represented by arrows. The 
constraints on the kinetics during each stage are indicated beside the arrow, (i) Slow nucleation 
stage (time scale r nuc ), in which critical nuclei of size m nuc form in a monomeric solution; this stage 
is found to be strongly hindered by kinetic barriers (dashed arrow), (ii) Fast growth stage, in which 
the nuclei grow from ?7i. nuc to an intermediate size m without changing their concentration. The 
growth may be either diffusion-limited (time scale r<j; —3/2 powerdaw relaxation) or kinetically 
limited (time scale t^; exponential relaxation), (iii) Bulk diffusion from the reservoir until the 
aggregates reach their equilibrium size m*; this stage is found to have a minor contribution to 
the micellization. (iv) Final relaxation of aggregate concentration through nucleation of additional 
micelles (also kinetically hindered; dashed arrow). 



The consequent prediction is that the formation of micelles in such open monomeric 
solutions may be suppressed for a macroscopic time. In fact, this behavior is regularly 
manifest in applications involving micelle-enhanced ultrafiltration [39( and has been observed 
in dialysis experiments [40|, where the time scale of micelle formation was estimated to be 
1-10 hours. In both the ultrafiltration techniques and the dialysis experiment a micellar 
solution is forced through a membrane, whose pores are smaller than the micelles. The 
surfactant solution on the other side of the membrane remains monomeric for a macroscopic 
time despite its contact with a micellar solution above the cmc. To our best knowledge the 
analysis presented above provides the first quantitative account of this regularly observed 
behavior. 

Apart from the aforementioned strong assumption of time-scale separation, the main 
shortcoming of our model is its mean-field character. We have assumed that the kinetics in 
the surfactant solution can be described within a representative subvolume, V±, containing a 
single aggregate and being uncorrelated with the other subvolumes. Upon closer inspection, 
in fact, we find that V\ for a closed system typically contains ~ 10-10 2 surfactant molecules, 
which is comparable to the aggregation number. Hence, correlations among such subvolumes 
are to be expected as the micelles nucleate and grow. Another important mean-field aspect 
is our description of the state of the system as a deterministic point on the free-energy 
landscape, and its kinetics — as a sharply defined path on that landscape. In practice, and 
particularly close to the cmc, the system should be more accurately described by stochastic 
distributions, with polydispersity and occupancies of both the monomeric and aggregated 
states [36| . Nonetheless, we do not expect these approximations to qualitative change the 
main results presented here. 
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Appendix 

In this appendix we calculate the asymptotic time dependence of the micellar size, m(t), 
in a diffusion-limited growth. The equations to be handled are eqs [TBH2"U1 

To leading order at long times we can substitute in eq [20] m(t) ~ fh, turning the 
boundary condition far away from the micelle into $i(r — > oo,t) = $i. We now define 
ip(r,t) = $i(r, t) — $i, so that ip(r — > oo,t) = 0, and introduce Laplace-transformed vari- 
ables, ?p(r,s) = J °° dte~ st ip(r,t), m(s) = J °° dte~ st m(t). The diffusion equation, eq CEE] is 
then rewritten as 




# = D-— r 2 ^f , (Al) 

where we have assumed ip(r, 0) = 0, as the accurate initial profile should not affect the 
long-time asymptotics. The boundary conditions, eqs [19] and [201 transform to 



AnR 2 dip 

sm - m nuc = D — — 

ncr dr 



(A2) 

r=R 

^(r ->• oo,t) = 0. (A3) 
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The solution of eqs IA1HA3I is 

na 3 srh — m nuc e ~( s / D ) 1/2 ( r ~ R ) 
AtiD 1 + R(s/Dy/ 

from which we get 



W S ) = - TIHi , o/./^M/ 2 ' ( A4 ) 



^ s) = - AnDRl + R(s/Dy^ (A5) 

T he lim it t — > oo corresponds to s — )■ 0, at which srh — m nuc ~ rh — m nuc = Am. Inverting 
eq lA5l back to real time and taking the limit t — > oo, we find 

na 3 Am , , , 

^•^^"-spSF 5 - (A6) 

which yields eqEUfor Td. 
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